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Charged Higgs boson may be produced at hadron colliders via the decay of top quarks. At the 
LHC, The large production cross section of top quarks pairs will make it one of the earliest channels 
allowing to search for physics beyond the standard model. Assuming the decay H + -¥ t + u, previous 
searches have so far considered only hadronic r decays. Here we show that by using appropriate 
kinematical variables (transverse mass observables), leptonic r decays can be used as well, for both 
observing and possibly measuring the mass of the charged Higgs boson. This can increase the overall 
experimental sensitivity to charged Higgs bosons at hadron colliders. 
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I. INTRODUCTION 

Charged Higgs bosons arise in models with extended 
Higgs sector such as two Higgs doublets models (2HDM) 
and in particular the Minimal Supersymmetric Standard 
Model (MSSM). If the charged Higgs mass is below the 
top mass, It is expected to be produced at hadron col- 
liders dominantly via top quark decay, due to the large 
Yukawa coupling of the top. At the LHC, the large pro- 
duction cross section of top quark pairs means that the 
existence of charged Higgs bosons could be probed with 
early data, with a potential of being among the first ev- 
idences of physics beyond the standard model. Here, we 
present methods for observing such a charged Higgs and 
possibly measuring it's mass, assuming that it decays to 
a r lepton and a neutrino, H + — > t + v, and using the 
leptonic decay modes of the r, either r + — > e + v e v T or 
t + — > /i + ^z7 T . The decay mode H + — > t + v dominates 
in the MSSM for tan/3 > 1 and is therefore a favorable 
search channel. Previous searches at the Tevatron [HQ 
considered only hadronic r decays, since those allow the 
t lepton to be identified and an excess of r leptons over 
the standard model prediction can provide a direct evi- 
dence for a charged Higgs boson. On the other hand, the 
leptonic decay modes involving an isolated electron or 
muon have experimental advantages with regard to trig- 
gering the event and suppressing backgrounds, and they 
are expected to be associated with less systematic un- 
certainties compared to the hadronic mode. This makes 
these channels particularly important for early stages of 
the LHC. 



The event topology we consider here belongs to a wide 
class of events where the presence of undetectable parti- 
cle^) in the final states prevents full reconstruction of the 
sought after particle's mass. It is often possible in such 
cases to define kinematical variables that are bounded 
by the unknown mass, and therefore it's value can be 



inferred from an edge in the distribution of that vari- 
able over many events. We generally refer to such vari- 
ables as transverse mass observables, following the most 
known example which is the W transverse mass in the 
process W — > Iv. Many generalizations to more com- 
plicated processes have been studied inpast years, such 
as the ttit2 variable and others [4] [5] 6]. In the fol- 
lowing sections we apply this approach to charged Higgs 
searches in tt events. We consider both semi-leptonic and 
di-leptonic final states, and show that such observables 
could be constructed in those cases and provide valuable 
information for observing and measuring the mass of the 
charged Higgs boson. 



Event simulation 

To study the performance of the proposed variables we 
use MADGRAPH @] for event generation, interfaced to 
PYTHIA [8( for showering and hadronization. t decays 
are handled by TAUOLA Detector effects and jet 
clustering are done with the PCS [l(| fast simulation 
package , configured with LHC parameters. All samples 
are generated for a center of mass energy of 14 TeV. 



II. SEMI-LEPTONIC tt EVENTS 

In this section we consider events in which one of the 
top quarks decays to a charged Higgs (t — > bH) while the 
second top quark decays hadronically and is therefore 
fully reconstructable (t — > bW — > bjj). This is the most 
useful channel for top reconstruction since the missing 
energy and isolated lepton together with b-tagged jets 
provide a good rejection against backgrounds. The main 
background to this process is the standard model decay 
t — >• Wb with a subsequent leptonic decay of the W, either 
directly W — >■ £v or via a r : W — > tv — > ivvv. It is 
natural for this channel to reconstruct the W transverse 
mass, defined as 11111 
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where p^.,p™ lss are the lepton and missing transverse mo- 
menta respectively and 4>£^ m iss is the azimuthal angle be- 
tween them. For direct W decays this transverse mass 
has a Jacobian peak near the W mass, hence it can pro- 
vide separation between the two W decay modes, as well 
as charged Higgs decays (see Figure [T]) . 
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FIG. 1: The W transverse mass distributions for leptonic W 
decays and charged Higgs boson decays with m H + = 130, 150 
GeV. all distributions are normalized to unit area. A mini- 
mum cut of 20 GeV for the lepton transverse momentum is 
applied. 

In principle, the presence of charged Higgs boson de- 
cays will be manifested as an excess of events in the lower 
region of the transverse mass distribution. In practice 
however, to observe such an excess would require an ac- 
curate modeling of the shape of the distribution which 
depends on the missing momentum resolution. This may 
be hard to achieve. A second drawback of this approach is 
that almost no information about the mass of the charged 
Higgs boson can be obtained. Still, by requiring rn^ to 
be sufficiently below the peak region, one can reject a 
large part of the background and enhance the signal con- 
tent of the sample. 

It is useful to define the W transverse mass (QJ equiv- 
alently as a minimization result, i.e. 

(m^) 2 = min [(p e + p mlss ) 2 } (2) 

{or«"?=o } 

where p ^p miss are the lepton and neutrino (missing) four 
momenta, respectively. It is straightforward to show that 
the above requirement indeed leads to the expression (|TJ) 
(see a similar derivation in appendix A). The definition 
(0) makes explicit the fact that m^ < mw, he- "m^ is 
bounded by the W mass. Furthermore, it implies that 
m^ is the best bound that can be placed on mw> per 
event. 

For the case of a leptonic t decay (either from a W 
or charged Higgs boson) the constraint used in ([2]) is not 
valid since the missing momentum comes from three neu- 
trinos and therefore (p miss ) 2 ^ 0. However, if one of the 



two b-jets in the event could be associated with the lep- 
tonically decaying top quark then the on-shell constraint 
for the top quark could be used, and one would define 
analogously a "Charged Higgs transverse mass" in the 
following way: 

(mf) 2 = max [(p e + p mlss ) 2 } (3) 

note that in this case the transverse mass is defined by 
maximization of the invariant mass, since it is bounded 
from above by the top quark mass. The charged Higgs 
transverse mass by definition therefore satisfies m H + < 
< m top , where is the true charged Higgs mass 

and m t op is the nominal value of the top mass used in the 
constraint. The explicit expression for it can be easily 
derived and is given by: (see appendix A) 

(m? f = ^ml p + Wt + P% + PT SS ) 2 ~ Pt) " (pt + Pt"') * 

(4) 

It should be noted that top quark decay can involve 
gluon radiation 12]. In such case the constraint used in 
is not exact. If the top quark emits a gluon and turns 
off-shell before it's decay, then the invariant mass of the 
decay products is smaller then the top pole mass 

{p e + p b + p miss ) 2 = m* top < m top (5) 

however, the charged Higgs transverse mass |0J| is a 
monotonically increasing function of the input value of 
mtop and therefore the following relation still holds: 

{mtop) > m" (m* top ) > m H+ (6) 

i.e., mif retains the property of being bounded by the 
charged Higgs mass also in the presence of gluon radia- 
tion. The values of however will be shifted upwards, 
making the Jacobian peak less pronounced. 



Event reconstruction 

To reconstruct the charged Higgs transverse mass, one 
must correctly associate a b-jet to the leptonic top decay. 
For the purpose of our study we consider events where 
exactly two of the jets are b-tagged, such that there are 
two possible choices. The selection cuts we apply are the 
following: 

• At least four jets with pt > 20 GeV and \r)\ < 2.5 

• Two of the above jets are b-tagged 

• Exactly one isolated lepton (e or /j,) with pr > 20 
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To choose between the two possible b-jet assignments, we 
first consider the invariant mass mf b = (p +p b ) 2 . For 
the correct b-jet, this is kinematically bounded by 

m 2 a < m 2 top - m 2 w (7) 

note that this bound holds for both W and charged Higgs 
boson decays, since we are assuming m#+ > mw [131 • In 
case that one of the lepton-b pairings violates the above 
condition, the opposite pairing in selected. In events 
where both pairings satisfy ([?]). we further consider the 
top quark transverse mass : 



p*$-p% iss ) 



(8) 



which is required to satisfy m T < m tov . Again, if one of 
the pairings violates this condition then the opposite one 
is selected. Finally, if this is not the case, the hadronic 
top is reconstructed using two of the light (non b-tagged) 
jets plus one of the b-jets. The b-jet that results in the 
invariant mass rrijjb being closest to the nominal value 
of rritop is associated to the hadronic top, and thus the 
other b-jet is associated to the leptonic top. 

Figure ^ shows the charged Higgs transverse mass 
distributions after selection and reconstruction. The 
power of discriminant based on the mass dif- 

ference is clearly observed. 
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FIG. 2: The charged Higgs transverse mass distributions for 
leotonic W decays and a charged Higgs decays with m H + = 
130, 150 GeV. all distributions are normalized to unit area. 

A useful property of is that the background dis- 
tribution is independent of the W decay mode (direct or 
r- mediated). This is easily understood from (|2J), which 
depends on the lepton and missing transverse momenta 
only via their sum pi, + p!£ lss = p^ . 

In figure |3] an example pseudo-experiment containing 
both the main background and a charged Higgs sig- 
nal with mjj+ — 130 GeV is shown. The number of 
events correspond to an integrated luminosity of 1 fb _1 
at y/s = 14 TeV, and the branching ratio Br(t — > H + b) 
is assumed to be 15%, roughly the current upper limit set 
by the Tevatron searches [2]. The statistical significance 
of discovery in this case exceeds 5cr, demonstrating the 
potential power of this variable. 
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FIG. 3: Example pseudo-experiments (points with error bars) 
on top of the signal and background distributions, assuming 
Br(t — > H + b) — 15% and m H + = 130 GeV, for an integrated 
luminosity of 1 fb _1 . 



III. DILEPTONIC tt EVENTS 



In this section we generalize the above procedure for 
dileptonic tt events, as depicted in Figure |U In this case 
the final state includes two leptons and missing energy 
on both sides of the event, making the reconstruction of 
the event more complicated. 




FIG. 4: A Feynman diagram illustrating the event topology 
of a di-leptonic tt event involving a charged Higgs decay. 



To define an analogue to the charged Higgs transverse 
mass of the previous section, we follow the same pro- 
cedure, i.e. maximize the invariant mass subject to all 
available constraints that are given by the known W and 
top quark masses, as well the conservation of momentum 
in the transverse plane. We have the following set of six 
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constraints, corresponding to the notation of Figure [U : 
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(9) 



Here p ff and p Vl represent the unknown quantities of 
the event. Note that we have not specified constraints for 
the particles that are the decay products of the charged 
Higgs, since this would not supply us with additional in- 
formation that can be use to constrain the charged Higgs 
mass. The system of constraints © therefore amounts to 
two free parameters (compared to only one in the semi- 
leptonic case) over which we maximize the charged Higgs 
mass, to obtain the variable m^ 2 



of shown in Figure O The effect of having addi- 

tional free parameter clearly worsens the situation com- 
pared to the semi-leptonic case, however there is still a 
clear discrimination between the W and charged Higgs 
bosons. 
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If we choose one of the free parameters to be the z- 
component of the charged Higgs momentum, we can im- 
mediately perform the maximization over it using the 
result of the previous section, and so m^ 2 is equivalent 
to 



l T2 



max 

{constraints} 



4{p" + )\ (11) 



with 



(mr(pf + 



The maximization over the remaining parameter needs 
to be performed numerically. A computation procedure 
which allows for this maximization to be easily performed 
is given in appendix B. 

Selection and reconstruction 

The selection cuts we apply in this case are similar 
to those used for the semi-leptonic channel, except that 
we require an additional isolated lepton with pt > 15 
GeV. The two leptons and two b-jets need to be paired 
according to the two corresponding top decays in order 
for mtf 2 to be calculated. To do this we calculate the 
invariant mass rnj b for each possible pair and require it 
to satisfy the bound ©. We keep only events where this 
requirement resolves the ambiguity, i.e. events where ex- 
actly one pairing satisfies (O . This keeps only about 50% 
of the events, however the accuracy of the assignment is 
high, at about 95%. The calculation of m^ 2 further re- 
quires assigning one of the two leptons to the charged 
Higgs decay. Here we always choose it as the one with 
lower transverse momentum. The resulting distributions 



FIG. 5: distribution of for dileptonic tt events containing 
a charged Higgs with m H + = 130, 150 GeV (red and green) 
and background events with a W boson decays (black) 



In figure [S] an example pseudo-experiment containing 
both the main background and a charged Higgs signal 
with m H + = 140 GeV is shown. Here the number of 
events correspond to an integrated luminosity of 10 fb _1 
at y/s = 14 TeV, and the branching ratio Br(t — > H + b) 
is assumed to be 15%. In this case, since the kinematic 
edge is less pronounced, to observe a charged Higgs using 
the m^ 2 variable requires a rather accurate knowledge of 
the background shape. The excess of events above the 
charged Higgs mass would then provide an indication for 
the charged Higgs in this channel, complementary to the 
semi-leptonic one. 
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FIG. 6: Example pseudo-experiments (points with error bars) 
on top of the signal and background distributions, assuming 
Br(t -s> H + b) = 15% and m H + = 140 GeV, for an integrated 
luminosity of 10 fb . 
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IV. CONCLUSIONS 

We have demonstrated how the leptonic decay chan- 
nel of the t could be used as a search channel for a 
charged Higgs boson at hadron colliders. The charged 
Higgs transverse mass defined for semi- leptonic ti events, 
as well as the analogous generalization for di-leptonic 
events, are sensitive to the charged Higgs mass and could 
therefore be used as discriminating variables that would 
give evidence to such a particle, and a way to measure 
it's mass. This would provide an important complemen- 
tary measurement to methods based on hadronic r iden- 
tification, and an enhancement to the overall discovery 
sensitivity. 

Acknowledgements 



plugging this into Eq. (A. 4) we obtain for A : 

1 - A = ^ (A7) 

where we have approximated the b-quark to be massless. 
Plugging (A. 7) and (A. 6) into (A. 3) we obtain the final 
result: 

(m¥f = {y/m* top + (f T + p% + P ™) 2 -Pt) 2 -(/t+pT SS ) 2 

(A8) 

It should be noted that exactly the same procedure could 
be applied to the case of a leptonic W decay W — > Iv , by 
replacing the constraint (A. 2) with (p™ ss ) 2 = . This 
will result in the usual expression for the W transverse 
mass. 
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Appendix A: Derivation of ra? 



Appendix B: Computation procedure for m T2 

In this section we describe the computational proce- 
dure by which the maximization of Eq. [TT]is performed. 
The transverse component of the charged Higgs momen- 
tum ptf is constrained by the set of equations 



In this appendix, we derive the expression for the 
charged Higgs transverse mass ([%])■ Our goal is to maxi- 
mize the expression of invariant mass: 



m' H+ = (p e +p miss f 

with respect to the two unknown components of p r ' 
while holding the constraint: 



(Al) 



l top 



(p n 



b-,2 



(A2) 



For convenience, we adopt the notation p« = (E,p z ) , 
satisfying p 2 = E 2 — p 2 z ,and pr = (p x ,p y ) with p\ = 
p 2 +p 2 ■ Equation (A.l) and (A. 2) can be re-written as: 



in 



h + = ip\ + P" liss r - (fx + p? ls T (A3) 



m\ op = (pjj + P 'f ss + P \f- (f T + Pt* ss + P%? (A4) 

To maximize Eq. (A. 3), we introduce a Lagrange 
multiplierA and differentiate with respect to p^ lss : 







dp™ ls: 
This gives 



{{p\ + p^ ss Y - (f T + i% lss Y) = (A5) 



Pt + -pS + +Pt° =Pt hiss 



{p vi y = o 

(p £ - +p^f = m 2 w 



(Bl) 
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We wish to obtain a parametrization of p^ + as a func- 
tion of a single variable such that the above constraints 
are satisfied. To do this we rewrite the last two equations 
as the following linear matrix equation: 



(B2) 



where 

Ax 

and 



E -p z 



{%). 



-P* ~Py 

i~ t~ 

-Px -Py 



(B3) 



The zero-mass constraint of the neutrino is given in 
this notation by 



where g is the 1+1 Lorentz metric, g = (J 
using (|B2|) this can be written as: 



(B4) 



p miss = __ p t i _ p b 



(A6) 



C T (1- A T gA)t + 2£ T A T gfh-rh T gm = (B5) 
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where A = Ai _1 A2 and in = Ai~ m. 

The matrix 1 — A T gA is symmetric and can be shown 
to be positive-definite. It can therefore be decomposed 
in the following way: 

1 - A T gA = L T L (B6) 

e.g. via a Cholesky decomposition. Using this we ob- 
tain 

|L£ +u\ = \Ju T u + m T gm = y/Q (B7) 

where lu = (L T ) _1 A T grh. 
This is solved by 

e = L- 1 ( v ^f(0)-c) (B8) 



where f is the unit vector in the transverse plane, 

We can therefore express as a function of <j): 

Pt + = Pt" ss + PS + ~ m (B9) 

such that 

m" 2 = max [rrtr {4>)\ (BIO) 

The maximization over <j) can now be performed by 
any standard function minimization algorithm, or by 
scanning over the range < < 2ir. The result is the 
generalized charged Higgs transverse mass. 
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